Quantum states experimentally achieving high-fidelity transmission over a spin chain 
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A uniformly coupled double quantum Hamiltonian for a spin chain has recently been implemented 
experimentally. We propose a method for the determination of initial quantum states that will 
provide perfect or near-perfect state transmission for an arbitrary Hamiltonian including this one. 
By calculating the eigenvalues and eigenvectors of a unitary operator obtained from the free evolution 
plus an exchange operator, we find that the double quantum Hamiltonian spin chain will support 
a three-spin initial encoding that will transfer along the chain with remarkably high fidelity. The 
fidelity is also found to decrease very slowly with increasing chain length. In addition, we are able 
to explain previous results showing exceptional transfer using this method. 
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Introduction. — Quantum information processing 
(QIP) often requires the transfer of known or unknown 
quantum states from one subspace to another within 
an information processing device. In recent years, the 
quantum spin chain has become a prime candidate 
for quantum communication purposes such as these 
[ll-Q. In the simplest configuration, where the nearest 
neighbor couplings are considered to be equal, perfect 
state transmission is typically not possible between two 
single spin processors within a linear chain. In other 
words, there is typically a non- vanishing probability that 
the initial excitation amplitude can be found outside 
the receiving spin location Q at any given time. In 
principle, however, perfect state transfer (PST) can be 
realized by properly engineering the couplings between 
neighboring sites ^Sj]. High fidelity state transmissions 
can also be obtained using weakly coupled external 
qubits 0] , modifying only one or two coup lings [j , Q , 



IC 



or by encoding the states using multiple spins 
[13I [3] , a class of states were found to transfer very well 
across long XY coupled spin chains. The existence of 
PST has also been established for a variety of interacting 
media, including, but not limited to, the spin chain 
model [l5| . Recently, exact state swap through a spin 
ring has been investigated. It was shown that there is a 
straightforward approach to calculating the probability 
of the occurrence of an exact state swap |16| . 



The schemes developed in Ref. [15| prompted the fol- 
lowing question. Given an arbitrary spin chain Hamil- 
tonian, can we find initial states which can be used to 
enable high-fidelity state transmission? In this letter, we 
answer this question and show that for a uniformly cou- 
pled chain, there exists a particular state which reliably 
transfers quantum information over large distances. We 
use a multi-spin encoding scheme and find the existence 



of a three-spin encoding which can provide reliable state 
transmission. In this case, the encoding and decoding 
processes can also be realized easily [13|. This report 
is therefore important from an experimental perspective 
due to the ease of implementation which is typically fa- 
vorable. 

The method for identifying high-fidelity states. — Con- 
sider a spin chain consisting of N sites which evolves 
according to some Hamiltonian _ff in a single excitation 
subspace. Suppose for the moment the initial state of our 
system is |*(0)) = |1) = (g) |101...1> (g) \0)g, where 
A and B denote separate processors. After the system 
evolves, the state at time t will be 



\^{t)) = U{t)\l) ^ exp{~iHt)\'if{0)) , 



(1) 



where h is taken to be 1 throughout. Suppose that at 
some time r PST occurs, then 



\^{T))=U{r)\l) = \N) 



(2) 



where |N) — |0)^(8) |101...1)(g) We can use a permu- 
tation operator Pab to swap all states in A and B, then 
the quantum information can be transferred from A to 
B. The permutation operator can be expressed as: 



^^s=L A\l3A){aA\^\aB){(3B\), 

^ — ^af3 



(3) 



where a,/3 = 1,2,..., 2'^ represent the standard basis 
for the k qubits located in processors A and B. Clearly 
P\g = Pab and P^^ = 1. \a{l3)A(B)) refers to a state 
|a(/3)) in processor A (B). From Eq. ([2]) 



Then 



UiT)\l)^PAB\l) 



PabU{t)\1) = W{t)\1) = \1) 



(4) 



(5) 



2 



We introduce the unitary operator W{t) = PabU{t). 
From Eq. ([5|), if the state |1) is an eigenvector of the 
operator W at time r, PST occurs. The eigenvectors of 
W reveal information about the possibihties of a specific 
state transmission. The problem of solving Schrodinger's 
equation now becomes a standard eigen-problem of the 
operator W. 

Since W(t) is a unitary operator it has a complete set 
of orthonormal eigenvectors {|\E'm(0))}T corresponding to 
eigenvalues {Em}T, 



Wir) |*™(0)) = Em |*™(0)) . 
This can also be written as 

C/(t)|*,„(0)) = S„P)sl*„,(0)) 



(6) 



(7) 



where C/(r) |^'m(0)) is the wave function [^^(t)) of the 
system which was initially prepared in the eigenstate 
|*„(0)). If |*m(0)) is a product state 



|*™(0)> = \A) (g> \C) 



(8) 



with \A) describing the state of processor A and |C) de- 
scribing the rest of the system, we can then obtain 



For the single excitation subspace, if one of the eigen- 
vectors l^'m(O)) = |1) at time r, PST occurs. If the 
eigenvectors are degenerate, an arbitrary linear superpo- 
sition of these degenerate states is also suitable for PST. 
Suppose there are L degenerate eigenvectors |^'i(0)) {I — 
1, 2, ...L), which have common eigenvalues E^. The state 



(10) 



1=1 



is an eigenvector of W{t), where Ci is an arbitrary num- 
ber. Our analysis describes a method for finding a state 
which can realize PST. (Note that these states are not all 
unique.) For a given Hamiltonian, if we initially prepare 
the state |5'(0)) as an eigenvector of the operator W{t), 
then after time r PST occurs. 



FIG. 1: Schematic of our quantum transmission protocol: (a) 
single-site encoding (b) three-site encoding. 



For state transmission from one end spin 1 to an- 
other end spin N, the exchange operator is given by 
Pi = |1) {N\ + \N) (1| + Po which is shown in Fig. [fl 



where Pq = J2j \j) 010 Ij-^)- We now considered a 
particular Hamiltonian, but emphasize that our method 
can be used for any Hamiltonian not only this example. 

An experimentally implementable Hamiltonian - Now 
consider the recently implemented Hamiltonian called a 
double quantum (DQ) Hamiltonian [17|: 



N-l 



(11) 



where Ji,i+i denotes the coupling between sites i and 
i + 1. This nearest neighbor coupled one-dimensional spin 
chain can be experimentally implemented using solid- 



|*™(r)) - E„,P\s \A) ® \C) = Em \B) ® \C') . (9) tity, [(E 



state nuclear magnetic resonance [17H19| in ^^F spins in 
a crystal of fluorapatite ((FAp-Ca5(POa)3)F) ^ HI]. 
The system described by Eq. (ITT|) will exhibit free evo- 
lution such that the evolution operator at time r will 
be U{t) =exp[— iriJ]. We can diagonalize the Hamil- 
tonian H such that Hd = W'^HW in the single exci- 
tation subspace. The evolution operator can therefore 
be expressed by U{t) = Wex-p[~iTHd]W and the N 
eigenvectors of W{t) can be obtained as a function of 
T. Furthermore, we consider a natural configuration for 
a DQ Hamiltonian with open ends. The z-component 
of the total for the staggered spins is a conserved quan- 



iGodd * 
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Zi),H] = 0. For simplicity, 
we will only consider the single excitation subspace of the 
full Hilbert space. In this case the total number of flipped 
spins is one. The basis for this subspace will be denoted 
as |j) which indicates that, after flipping, the even (odd) 
site spins all of the spins reside in the |0) (|1)) state except 
for the spin at site j which is in the |1) ( |0))state. For 
example, in a = 5 site chain, the single excitation sub- 
space wiU be spanned by |1) = |11010) , |2) = |GGG10),..., 
etc. If we flip the even numbered states we find that the 
total up spin is actually one. We will use this description 
throughout this paper. 

Example I: nonuniform couplings — We will consider 
several different coupling configurations with the poten- 
tial for high-fidelity state transmission and the best re- 
sults will be provided at the end of our analysis. First 
as an example, we consider two pre-engineered couplings: 
(1) weak couplings at both ends, where Ji_2 = Jn-i.n — 
Jo and Ji,i+i = J elsewhere. (2) couplings termed PST, 
where Ji,i+i = \/i{N — i). It is already known that high 
fidelity (Pmaa ~ 1) state transmission for the first con- 
figuration ]^ and perfect fidelity {Fmax = 1) for the 
second configuration can be gained in a spin system 
Here we will use these two kinds couplings to show the 
applicability of our methods. 

For a five-spin system with weak couplings at both 
ends, we take Ji^2 = '^4,5 — O.IJ. J equals -1 elsewhere. 
The eigenvalues and eigenvectors of the operator W{t) at 
an arbitrary time r can be obtained numerically. We will 
consider those which span the single-excitation subspace. 
In Table m we plot the results for r = 31. The first col- 
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|1) 12) 


|3) 


|4) 


|5) 


1 (0.999,0.03) 0.035 -0.500 


0.705 


-0.500 


0.035 


2 (0.999,-0.03) 0.035 0.500 


0.705 


0.500 


0.035 


3 (1.000,0) 0.398 


-0.066 





0.914 


4 (1.000,0) 0.999 


-0.050 





-0.003 


5 (-1.000,0) 0.707 





-0.707 






where 



TABLE I: The complex eigenvalues and corresponding eigen- 
vectors of the operator W{t) at r = 31 in a spin chain where 
the weak coupling conditions Ji,2 = J4,5 ~ O.IJ are satisfied. 
We take N=5, J = -1. 



umn labels the complex eigenvalues while the remaining 
columns are associated with the amplitudes of the states 
at the top of each column. The coefficients in columns 
2-6 could be complex numbers but our results show that 
the imaginary part always equals zero, so we take them 
to be real throughout. The same meaning holds for Ta- 
bles m ini IIVI and |Vl For example, at the line labeled 
with a 4 the eigenvalue is (1.000,0) and the eigenvector 
is 0.999 |1) - 0.050 |3) - 0.003 |5). The state |1) closely 
approximates this eigenvector and it can be written as 
the aforementioned product state. If we use the state 
|1) as the initial state of the whole system, then at time 
T = 31, the system will closely approximately the state 
IN). 

We use the fidelity between the received state and the 
ideally transfered state, F = ^/{^{0)\p{t) |$(0)) as a 
measure of the quality of the transfer. Here |$(0)) is 
a state at the receiving end which has the same form as 
the state initially prepared by the sender. p{t) is the re- 
duced density matrix of the receiver's spin at time t and 
is obtained by tracing over all but the receiver's sites. In 
Fig. [2] (a) we plot the fidelity versus time t for the weakly 
coupled chain. The initial state is |1) has the maximum 
fidelity, = 1 at time t = 31. 



Il) 


|2) 


13) 


|4) 


1 (0,1.000) 0.707 








-0.707 


2 (0.018,-1.000) 


0.707 


-0.707 





3 (-0.009,1.000) 0.612 


-0.354 


-0.354 


0.612 


4 (0.028,1.000) 0.354 


0.612 


0.612 


0.354 



TABLE IL The complex eigenvalues and corresponding eigen- 
vectors of the operator W{t) at t — 3.14 in a A'^=4 site spin 
chain with couplings given by Ji,,;+i = \/i{N — i). 



Next we discuss PST. For the simple case of = 4, 
the results at time r = 3.14 are listed in Table ITTl None of 
the eigenvectors can be written in the form of a product 
state |4'„(0)) = 1^) (g) |C), but the eigenvalues of 1, 3, 4 
are roughly degenerate. Consider the superposition 



1*1(0)) = 
1*3(0)) = 

1*4(0)) = 



(|1)-|4)), 
(|1> + |4)) 
(|1> + |4)) 



(|2) + |3)), (13) 
(|2) + |3)). 



The state 1 1) at site 1 can be transferred exactly to site 
4 at time r = 3.14. In Fig. ^h) we plot the time evo- 
lution of the fidelity when transferring a state |1) from 
site 1 to 4. We also see that at time t = 3.14 the fi- 
delity is nearly 1. These examples illustrate the validity 
and practicality of our method while providing a general 
method to obtain the results. 





FIG. 2: The fidelity as a function of time for (a) N — 5, weak 
coupling conditions for Ji,2 = J4,,5 = O.IJ (b) A'^ = 4, channel 
coupling conditions Ji,i+i = \/i{N — i). See text for more 
details. 

Example II: uniform couplings. — Consider the most 
natural configuration for a spin chain; a uniformly cou- 
pled spin chain. We take the ferromagnetic coupling 
Ji.i+i = J = — 1. Note that PST is typically unattain- 
able in these systems using single-spin encodings [lo| . 
We will consider both single spin encodings as well as 
multi-spin encodings. For single-spin encodings, our cal- 
culations confirm that PST cannot occur in this model 
as is known [l6| . In table IIIII we list the the maximal Pm 
and the corresponding r for a N — 7 uniform chain for 
different eigenvectors. Here Pm is defined as the overlap 
between the eigenvector |*m(0)) and the initial state |1), 
i.e.,p™ = |(*™(0)|1)|. 

l*,n) |*i) WT) W) 1*1) W) W) WTT 



31.1 
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31 



18.6 35.6 30.6 



8.9 



V2|«'i(0)) 



1*3(0)) + J- 1*4(0)) = |1), (12) 



Pm 0.7071 0.6295 0.7062 0.6402 0.7068 0.6842 0.7071 

TABLE III: The maximal p,„ and corresponding values of r. 
N—7, the maximum values are found in a time interval [5, 40]. 

Now we will examine multi-spin encoding schemes. 
As an example, we first consider a three-spin encoding. 
Specifically, suppose we wish to transfer a state of the 
form |*(0)) = [a |100) -t- /3 |010) -I- 7 |001))a «) |00...0) (g 
|000)^. As shown in Fig. [Ub), we intend to transfer the 
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state of the first three spins to the opposite end. The 
results for a = 6 site chain are given in Table HVl for 
time T — 4.0. The eigenvalues of 1 and 2 are roughly 
degenerate and the approximate relation 



|*i(0)) + 1^1/2(0)) =- |1) + |3) 



(14) 



can be written in the form of a product state (— |110) + 
|011))a(»|10...1). The state (- |110) + |011))/v^ is there- 
fore suitable for transmission. We have also checked the 
case where N — 7 and find that at time r = 28.8 the 
above states can be obtained again. 







|1> 


|2) 


|3> 


|4) 


|5) 


|6) 


1 


(0.117,-0.993) 


-0.493 


0.005 


0.500 


0.500 


0.005 


-0.500 


2 (-0.117,-0.993) 


-0.493 


-0.005 


0.500 


-0.500 


0.005 


0.500 


3 


(0.252,0.968) 


-0.275 


0.590 


-0.275 


-0.275 


0.590 


-0.275 


4 


(-0.252,0.968) 


0.275 


0.590 


0.275 


-0.275 


-0.590 


-0.275 


5 


(0.544,0.839) 


0.421 


0.379 


0.405 


0.405 


0.379 


0.421 


6 


(-0.54,0.839) 


-0.421 


0.379 


-0.405 


0.405 


-0.379 


0.421 



TABLE IV: The eigenvalues and corresponding eigenvectors 
of the operator W{t) at t — 4.0 using a 3 spin encoding. Here 
iV=6. 





|1) 


|2) 


|3) 


|4) 


|5) 


1 (0.999,-0.025) 


-0.263 


-0.263 


0.850 


-0.263 


-0.263 


2 (-0.999,0.041) 


0.500 


-0.500 


0.000 


-0.500 


0.500 


3 (-0.997,0.076) 


0.500 


-0.500 


0.000 


0.500 


-0.500 


4 (0.997,0.076) 


0.500 


0.500 


0.000 


-0.500 


-0.500 


5 (0.998,0.067) 


0.425 


0.425 


0.526 


0.425 


0.425 



TABLE V: The eigenvalues and corresponding eigenvectors of 
the operator W{t) at t — 47.2 using a 2-spin encoding. Here 
N=5. 
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FIG. 3: (Color online.) Length dependence of the maximum 
fidelity achievable Fmax and the associated arrival times Tmax 
for the state (a) (- |110) -h |011) )/^/2 and (b) (|11) - |00))/V^. 
The time is searched within the interval [0, 50]. 

Table |V] lists the results corresponding to a = 5 
chain for the two-spin encoding. At r = 47.2, the relation 
1*2(0)) -I- |^'3(0)) = |2) - |1) approximately holds which 
can also be written in the form of the product state in 
Eq. ©. 

We have found some states realizing high-fidelity state 
transmission, for small N. Now we will check to see if 



the state (- |11) -I- |00))/V2 and (- |110) + |011))/%/2 
can be transferred with high fidelity across chains of 
arbitrary length N. In Fig. |3l we plot the maximum 
fidelity Fmax and the associated arrival time Tmax as 
a function of chain length N . The analytic expression 
with eigenvalues Em = — 2 Jcos[7rm/(Af + 1)] and eigen- 
vectors |^',„(0)) = ^/'2/{N+ l)EjSin(gmj) |j) are used. 
For practical implementation of our protocol, the maxi- 
mum fidelity is found in the time [0,50]. For the two-spin 
encoding, the high fidelity associated with short chain 
lengths cannot be acheived with increasing chain length. 
Fmax quickly decreases with increasing TV. However, this 
robustness can be observed even for long chains using the 
three-spin encoding. The fidelity is exceptionally large 
for a relatively long chain. Therefore, using this state, a 
high-fidelity state transfer can be gained. Fmax ~ 0.96 
for = 6 at t = 4.0, Fmax = 1-00 for iV = 7 at t = 28.8 
which agrees with our previous analysis 



13, 14 . Note 



that we only consider two and three-spin encodings here. 
For encodings using more than three spins, we conjecture 
that for odd spin encodings some states can be found to 
possess high-fidelity transmission even over long chains. 
From Fig. |3] (b) we find that the arrival time Tmax typi- 
cally increases with increasing chain length TV except for 
some deviation with small values of N . We also find that 
the Tmax associated with the three-spin encoding is a lit- 
tle longer than in the two-spin encoding case for N > 24. 
This suggests that encodings using larger Hilbert spaces 
require longer waiting times for the maximum fidelity. 

Conclusions. -In conclusion, we have introduced a 
method to find states which can be transmitted through 
spin channels with high fidelity. The method can be eas- 
ily implemented numerically and can be applied to A^- 
site encodings, with A^ arbitrary. Using our method we 
have provided examples for the DQ Hamiltonian which 
exhibit uniform and nonuniform exchange couplings. For 
the uniform chain, a 3-spin encoding (- |110)-|-|011))/\/2 
was found to exhibit high fidelity state transmission. Us- 
ing a simple similarity transformation [l3|, our results 
can be extended to the standard Heisenberg XY model. 
In this case we have provided an explanation for the ap- 
pearence of the class of_initial states which were pre- 
viously discovered 



13, 14 



These states are exceptional 
due to the fact that they use simple encodings and trans- 
fer extremely well. Our work therefore provides a new 
method, new results, and an explanation of previously 
known important results. 

This work is supported by NSFC (Grant No. 
11005099), Fundamental Research Funds for the Central 
Universities (Grant No. 201013037), an Ikerbasque Foun- 
dation Startup, the Basque Government (grant IT472- 
10) the Spanish MICINN (Project No. FIS2009-12773- 
C02-02). MSB is supported by NSF 0545798. 



5 



* This submission was written by the author acting in 
his own independent capacity and not on behalf of UT- 
Battelle, LLC, or its affiliates or successors. 

[1] S. Bose, Phys. Rev. Lett. 91, 207901 (2003);Contem. 
Phys. 48, 13 (2007). 

[2] L.-A. Wu, A. Miranowicz, X.B. Wang, Y.-X. Liu, and 
F. Nori, Phys. Rev. A , 80 012332 (2009). L.A. Wu, 
M.S. Byrd, Z.D. Wang, and B. Shao, Phys. Rev. A 
82,052339(2010). 

[3] H. Yung, Phys. Rev. A 74, 030303 (R) (2006). 

[4] M. Wiesniak, Phys. Rev. A 78, 052334 (2008). 

[5] M. Christandl, N. Datta, A. Ekcrt, and A.J. Landahl, 
Phys. Rev. Lett. 92, 187902 (2004). 

[6] Antoni Wojcik ct al., Phys. Rev. A 72, 034303 (2005). 

[7] S. Oh, L.-A. Wu, Y.-P. Shim, J. Fei, M. Friesen, and 
X. Hu, Phys. Rev. A 84, 022330 (2011), and references 
therein. 

[8] T.J.G. Apollaro, L. Banchi, A. Cuccoli, R. Vaia, and P. 

Vcrrucchi, Phys. Rev. A 85, 052319 (2012). 
[9] L. Banchi, T. j.G. Apollaro, A. Cuccoh, R. Vaia, and P. 



Vcrrucchi, New J. Phys. 13, 123006 (2011). 
[10] T.J. Osborne and N. Linden, Phys. Rev. A 69, 

052315(2004). 
[11] H.L. Haselgrove, Phys. Rev. A 72, 062326(2005). 
[12] J. AUcock and N. Linden, Phys. Rev. Lett. 

102,110501(2009). 
[13] Z.-M. Wang, C.A. Bishop, M.S. Byrd, B. Shao and J. 

Zou, Phys. Rev. A 80, 022330 (2009). 
[14] C.A. Bishop, Y.-C. On, Z.-M. Wang, M.S. Byrd, Phys. 

Rev. A 81, 042313 (2010). 
[15] L.-A. Wu, Y.-x. Liu, and F. Nori, Phys. Rev. A 80, 

042315 (2009). 

[16] B.Q. Liu, L.A. Wu, B. Shao, and J. Zou, Phys. Rev. A 

85, 042328 (2012). 
[17] C. Ramanathan, P. Cappellaro, L. Viola, and D.G. Cory, 

New. J. Phys. 13, 103015 (2011). 
[18] E. Rufeil-Fiori, CM. Sanchez, F.Y. Oliva, H.M. 

Pastawski and P.R. Levstein, Phys. Rev. A 79, 032324 

(2009). 

[19] W. Zhang, P. Cappellaro, N. Antler. B. Pepper, D.G. 
Cory, V.V. Dobrovitski, C. Ramanathan and L. Viola, 
Phys. Rev. A 80, 052323 (2009). 



